A corresponding-state approach to quark-cluster matter 
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The state of super-dense matter is essential for us to understand the nature of pulsars, but the 
non-perturbative quantum chromodynamics (QCD) makes it very difficult for direct calculations 
of the state of cold matter at realistic baryon number densities inside compact stars. The strong 
coupling between quarks might render quarks grouped in clusters, and at high densities but low 
temperature, it is conjectured that the quark-cluster matter could be in a condensed phase as self- 
bound objects. Nevertheless, supposing that the quark-clusters could be analogized to inert gases, we 
apply here the corresponding-state approach to derive the equation of state of quark-cluster matter, 
as was demonstrated for nuclear and neutron-star matter in 1970s. According to the calculations 
presented, the quark-cluster stars, which are composed of quark-cluster matter, could then have high 
maximum mass that is consistent with observations and, in turn, further observations of pulsar mass 
would also put constraints to the properties of quark-cluster matter. Moreover, the melting heat 
during solid-liquid phase conversion and the related astrophysical consequences are also discussed. 

PACS numbers: 26.60.-c, 97.60.Gb, 21.65.Qr 



I. INTRODUCTION 

The state of matter above the nuclear matter density 
po is still far from certainty, whereas it is essential for 
us to explore the nature of compact stars. At average 
density higher than ~ 2p , the quark degrees of free- 
dom inside would not be negligible, and historically such 
compact stars are called quark stars [U-IH • Although cold 
quark matter is difficult to be created in laboratories or 
studied by direct QCD calculations, some efforts have 
been made to model quark matter and quark stars,from 
MIT bag model to color super-conductivity model [I| . In 
most of these models, quark matter is characterized by 
soft equation of state, because the asymptotic freedom 
of QCD tells us that as energy scale goes higher, the 
interaction between quarks becomes weaker. 

However, at realistic baryon densities of compact stars, 
p ~ 2 — lOpoj the energy scale is usually below 0.8 GeV, 
which is much lower than the scale where the asymptotic 
freedom could apply. In contrast, the non-perturbative 
effect should be significant, making quarks to couple 
strongly with each other. Quark-clustering is then con- 
jectured to occur in cold dense matter inside compact 
stars, by condensation of quarks in position space due to 
the strong coupling The realistic quark stars could 
then be actually "quark-cluster stars" , and solidification 
could be a natural result if the kinetic energy of quark- 
clusters is much lower than the interaction energy be- 
tween them. The idea of clustering quark matter could 
provide us a way to understand different manifestations 
of pulsars [f| . 

How can then one model the equation of sate of quark- 
cluster matter? Due to the lack of both theoretical and 
experimental evidence, the hypothetical quark-clusters 



in cold dense matter has not been confirmed, and it is 
also difficult for us to derive the properties of quark- 
cluster matter by direct QCD calculations. Nonetheless, 
an empirical way was employed and discussed seriously 
in a corresponding-state approach to deduce the prop- 
erties of nuclear matter in 1970s [e.g. 0, H|- To estab- 
lish a model which could be tested by observations, we 
adopt an empirical method here, by analogizing quark- 
clusters to inert gases and applying the corresponding- 
state approach too. A quark-cluster is assumed to be 
usually colorless, just like an inert atom is electric neu- 
tral. The interaction between inert gas atoms is the re- 
sult of residual electromagnetic force, and similarly the 
interaction of quark-clusters could be seen as the result 
of residual strong force, both of which should be charac- 
terized by the short-distance repulsion and long-distance 
attraction. In this paper, we assume that the interaction 
between quark-clusters could be described by the same 
form as that between inert gas atoms, i.e. the Lennard- 
Jones potential, only with different parameters indicat- 
ing stronger interaction and larger densities. A special 
kind of quark-cluster, so-called if-cluster, was studied 
extensively Q, and the interaction between .ff-clusters 
could be Lennard- Jones-like under the Yukawa potential 
with a and to coupling. In addition, nucleon could also 
be considered as special 3-quark clusters, the interaction 
between which is found to be Lennard- Jones-like by both 
experiment and modeling [Ioj |. 



In fact, quark matter in Lennard- Jones model has been 
studied, where the equation of state is derived by s um 
ming the interaction energy of all quark-clusters 
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The so-called corresponding state approach we demon- 
strated in this paper, however, is an empirical one, to 
derive properties of quark-cluster matter by just a com- 
parison to the experimental data of inert gases, based 
on the law of corresponding states. The law of corre- 
sponding states was first proposed by de Boer [l2j, who 
found that the properties of inert gases, such as pres- 
sure and density, could be written in a reduced form. 
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After reducing to dimensionlcss terms, the experimen- 
tal data of various inter gases can be fitted in smooth 
curves with a single quantum parameter. If the quark- 
cluster matter is assumed to be similar to inert gases, 
the corresponding-state approach can be applicable to 
study the state of quark-cluster matter, without knowing 
its exact structure. Although interaction between quark- 
clusters may not be perfectly described by Lennard- Jones 
potential, the corresponding properties of quark-cluster 
matter could be in the same ranges, so our approach is 
to some degrees reasonable. 

With the same form of interaction, we can derive the 
equation of state of quark- cluster matter from empirical 
data of inert gases, by the corresponding-state approach. 
The masses and radii of quark-cluster stars can then be 
derived and compared with observations. We find that 
the maximum mass of quark-cluster stars can be well 
above 2M Q . The melting heat is also discussed, and it is 
shown that the solidification of newly born quark-cluster 
stars might explain the plateau of 7-ray bursts. 

This paper is arranged as follows. A brief introduction 
of the law of corresponding states is given in §11. The 
equation of state of quark-cluster matter and the mass- 
radius curve of quark-cluster stars are derived in §111 us- 
ing the corresponding-state approach. The melting heat 
of solid quark-cluster stars and the related astrophysical 
consequences are discussed in §IV. We make conclusions 
and discussions in §V. 

II. THE LAW OF CORRESPONDING STATES 

The law of corresponding states, advocated by dc 
Boer [l2| , shows that the equation of state of substances 
with same form of interaction can be written in a re- 
duced and universal form. Consider a group of sub- 
stances with the following properties: (1) the total po- 
tential energy due to interaction can be written as a sum 
of identical expressions tp{r ik )^ each of which depends 
only on the distance between two particles i and fc; 
(2) ip(r) = ef(r/cr), where / is a function same for all 
substances, and e, a are characteristic energy and length 
for different species. The macroscopic quantities, such 
as pressure P, volume V and temperature T, can be ex- 
pressed in dimensionless terms: 

P* = Pa 3 /e (1) 
V* = V/(Na 3 ) (2) 
T* = kT/e (3) 

Another dimensionless parameter is 

A* = h/ (ay/me), (4) 

corresponding to the de Broglie wavelength, which is con- 
structed to measure the importance of quantum effects. 
It can be proved that the reduced equation of states ex- 
pressed in dimensionless quantities is a universal relation 

P* = f(V*,T*,A*), (5) 



which is the formulation of the law of corresponding 
states [T^| . 

Despite the so-called universal equation of states is just 
formally written as eq. ([S]) , a formula that is difficult to be 
derived theoretically for most cases, it could be used to 
obtain information on the equation of state of a substance 
which we are unfamiliar with. For determined V* and 
T*, P* depends on the value of A*, and the P* - A* 
curve can be drawn using experimental data of laboratory 
substances. If the curve is smooth enough, the value of 
P* for unfamiliar matter at such a state can be predicted 
provided its A* is known. 

For some substances described by Lennard- Jones 6-12 
potential 

de Boer had determined e and a of noble gases and some 
permanent gases (r = a is the distance where (f{r) = 0, 
and e is the depth of potential well) [I2J . Then the exper- 
imental data of P* , T* or V* for different substances turn 
out to be smooth functions of A* as corresponding states. 
In Fig[TJ experimental data of the volume Vq at zero tem- 
perature and zero pressure, reduced to Vq = Vo/{Na 3 ), 
are plotted with A* for some substances. A smooth curve 
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FIG. 1. Experimental data of reduced volume Vo* and A* 
at zero temperature and pressure for different inert gases are 
shown by dots [13], and the fitted curve of eq.© is shown by 
solid line. Four cases A, B, C, and D are denoted by crosses, 
which will be explained in §111 A. 

can be drawn by fitting all the points, which forms the 
bases of our prediction via corresponding states law, and 
the formula for the fitted curve is 

F * = 0.57 + 9.45 x 10~ 5 (A* + 6.35) 444 . (7) 

Considering the property of short-distance repulsion 
and long-distance attraction shown by Lennard-Jones po- 
tential, we assume that the interaction between quark- 
clusters can also be described by this form. The dis- 
tinctions between quark-cluster matter and ordinary sub- 
stances should be a much deeper potential well (larger 
e) and higher density (smaller a). With the same form 
of interaction as that of inert gas, we could apply the 
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law of corresponding states to derive the properties of 
quark-cluster matter. If we find the quantum parameter 
A* corresponding to quark-cluster matter, then V * and 
other properties that vary smoothly with A* can be de- 
termined by simply looking at the experimental curves 
of that property vs A* for inert gases. 



III. THE STATE OF QUARK CLUSTER 
MATTER 



A. Parameters 

To apply the law of corresponding states to quark- 
cluster matter, we must determine e, a and the mass m 
of each quark-cluster first, m depends on the number of 
quarks N q and the mass of each quark m in one cluster. 
We give each quark a constituent mass and assume mo is 
one-third of the nuclear mass. N q is left as a free param- 
eter in this paper, and we choose two cases coming from 
the following considerations. A dihyperon with quantum 
number of AA {H dibaryon) was predicted to be a stale 
state or resonance and recently Lattice QCD sim- 
ulations have found evidence for its existence [HI, EH- 
Besides, an 18-quark cluster, i.e. quark-a, being com- 
pletely symmetric in spin, color and flavor space, was 
also predicted to exist [16[ . Consequently, we set N q = 6 
and N q — 18 for our calculations. 

As no experimental attempt has been made to get the 
values of e and a, we try to constrain their values by 
the surface density p s of quark-cluster stars. The tem- 
perature of quark stars can be approximated to be zero, 
and the pressure also reaches zero at the surface of stars. 
Given the value of V * , we can calculate the surface den- 
sity p s (rest- mass density). It is obvious that p s can be 
written as 

Ps = N ■ N q m /V , (8) 
and comparing cq.© with cq.® we can get 

p. = N q m /(V *a 3 ). (9) 

For certain values of N q , e and a, we can calculate A* 
of quark cluster matter by eq.Q, and Vq can be found 
according to the fitted relation eq.© of Vq*-A* curve, 
then we may determine p s using eq.©. In Fig(2J pairs of 
£ and a that correspond to the same surface density p s 
are plotted respectively for N q = 6 and N q = 18, where 
values of p s are chosen to be once, twice and three times 
of nuclear matter density po- The lines of e and a giving 
the same A* with values 1, 2 and 3 are also drawn here 
for a further limit. 

The surface density of quark stars is assumed to be 
in the range 1 < p s /po < 3. Quark-clusters could con- 
densate to form solid state like classical particles, so the 
quantum effects may not be large for quark-cluster mat- 
ter, then A* should satisfy A* < 2. We select four points 
numbered A, B, C and D representatively to deduce the 



3.5 
3.0 
2.5 
2.0 
1.5 
1.0 













~ l~ 


=1 


















■= '<■< 






- Nq 18 




w 














& 










Ps 


'Po=2 










t 












Ps 'Po =1 ■ 
































P»'Po =2 - 














'D 










•< 




r\ V! 




Ps'Pcr 3 












\^ 





10 



100 1000 

e (MeV) 



FIG. 2. Contour lines of surface density p s and A* , with solid 
lines representing N q — 6 and the dashed lines representing 
N q = 18, including p B /p = 1,2,3 and A* = 1,2,3. Four 
cases A, B, C and D are denoted by stars. 



equation of state for quark-cluster matter and then the 
mass-radius relation of quark-cluster star. The values of 
e, a and the resulting p s . A* at points A to D are given 
in Table IIII Al They are also plotted in Fig. [1] corre- 
sponding to four different cases for quark-cluster matter, 
in which the equation of state will be calculated respec- 
tively. 



Table DTDS 





N q 


e (MeV) 


a (fm) 


A* 


Ps/po 


A 


18 


40 


2.5 


1.05 


1.87 


B 


18 


100 


2.3 


0.72 


2.72 


C 


6 


150 


1.5 


1.56 


2.47 


D 


6 


200 


2.0 


1.02 


1.23 



Our choice of parameters e and a comes from the fol- 
lowing consideration. The depth of potential well for nu- 
clear matter is about 100 MeV, so it may be reasonable 
that £ = O(100 MeV). The average inter-cluster distance 
d at the surface of quark-cluster stars is given by 



.3 x 0.16^/po), 



N„ 



(10) 



With Ps/pa = 2, we get d 
d = 2.66 fm for N n = 18 



1.84 fm for N q = 



6 and 

Then a = 0(1 fm) as it 
should have the same order of magnitude as d. It can be 
seen that the selected parameters are consistent with the 
above estimation. 



B. The equation of state 

Given e and a, we can deduce the state of quark-cluster 
matter by a corresponding-state approach, in the zero 
temperature case. If we know the experimental P* — A* 
curve at a certain V* and zero temperature, we can find 
the value of P* corresponding to A* of quark cluster. 
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FIG. 3. Experimental data of P* and A* at zero temperature 
when V* — 0.88 are shown by dots [T3|, and the fitted curve 
is almost a straight line (solid line). The cases A, B, C, and 
D are denoted by crosses. 
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FIG. 4. The internal energy E per cluster for four groups 
of parameters A(red line), B(black line), C(blue line) and 
D(cyan line). The range of density n is from surface den- 
sity n s to the highest central density where the quark-cluster 
stars reaches the maximum mass. 



According to eq.([T]l and the number density of quark- 
clusters n = l/(V*a 3 ), the reduced quantities P* and 
V* can be converted to P and n, then we can get the 
pressure at a certain number density of quark-clusters. 
Combining this with the relation between mass density 
p (rest-mass density plus interaction energy density) and 
number density n of quark-cluster matter, the equation 
of state can be derived. 

To draw the P* — A* curve at different V*, we need to 
know the relationship between P* and V* of some sub- 
stances at zero temperature. For the lack of new data, 
we just use the data provided by de Boer in his subse- 
quent article [13], where values of P* and A* were given 
corresponding to different values of V* for various in- 
ert gases. Taking V* = 0.88 for instance, the P* — A* 
curve are shown in Fig|3] The data points are almost 
in linear relation, which makes our interpolation reli- 
able. The value of A* at point A is about 1.05, then 
we find P* w 25 from the P* — A* curve. The corre- 
sponding P and n to the reduced quantities P* and V* 
are P = 1.0 x 10 35 dyn/cm 2 , u/uq = 2.7 (no is the num- 
ber density of nucleons in nuclear matter). Thus we get 
P(n = 2.7no) = 1.0 x 10 35 dyn/cm for quark-cluster 
matter in case A. By taking different values of V* , pres- 
sure P at different densities can be determined in case A. 
The same procedure is also applicable to the other three 
cases. 

For each set of parameters, what we get is just a set of 
points in P—n diagram and not an analytic equation, and 
then we perform the curve fitting to get an approximate 
formula. The P — n relations derived from curve fitting 
are 

P = (2.99 x 10 41 n 5 63 - 1.60 x 10 34 ) dyn/cm 2 (11) 
P = (1.99 x 10 41 n 5 64 - 7.63 x 10 34 ) dyn/cm 2 (12) 
P = (8.10 x 10 38 n 5 24 - 1.69 x 10 35 ) dyn/cm 2 (13) 
P = (6.69 x 10 40 n 5 63 - 1.63 x 10 35 ) dyn/cm 2 (14) 

for A, B, C and D respectively, where n is in units of 
clusters/fm 3 . Certainly it is better to deduce the equa- 



tion of state from a border range of densities, making the 
extrapolation to be more accurate. Nevertheless, lacking 
in experimental data of laboratory substances, we can 
only make such an approximation at this stage. It is 
worth mentioning that the approximation will not have 
much influence on the following calculations of the mass- 
radius curves. 

According to P = n 2 4^ , where E is the internal energy 
per cluster, we can get 

r" P(n) 

E(n)-E(n s )= / -\±dn, (15) 

where n s is the number density of quark-clusters on the 
surface of stars. We may determine the value of E(n s ) 
from a corresponding-state point of view, and then the 
relation between E and n can be derived by the above in- 
tegral. Similar to V *, Uq = Uq/Ne can be approximated 
as a smooth function of A*, where is the internal en- 
ergy at zero temperature and zero pressure. From the 
data of laboratory substances [l2[ , we derive a fitted for- 
mula for Uq, 

C/ * = -8.72 + 4.91A* - 0.71A* 2 , (16) 

and E{n s ) is thus 

E(n s ) = U /N = U*e. (17) 

As both P(n) and E(n s ) are known, it is able to calculate 
E(n) from cq. (fT5")) . The results are plotted in Fig01 for 
four groups of parameters A to D, and we can see that the 
internal energy can be comparable to rest-mass energy at 
some densities. 

The mass density p consists of rest-mass density and 
energy density, 

p = n{N q m + E/c 2 ) 1 (18) 

then the equation of state for quark-cluster matter can 
be derived by combining P — n relation and eq. (|18p , and 
we show the results in Figf5J for the four groups of pa- 
rameters A to D. 



FIG. 5. Equations of states for the same four groups of 
parameters as in Fig. [4] 



Mass-radius relation 



Considering perfect fluid case and the general relativ- 
ity, the hydrostatic equilibrium in spherically symmetry 
is described by Tolman-Oppcnheimer-Volkoff equation, 



dP 

~d7 



Gm(r)p (l + ^r)(l + 



m(r)c 2 i 



2Gm(r) 



(19) 
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FIG. 6. The mass-radius and mass-central density (rest-mass 
density) curves, and different parameters are distinguished by 
their colors as the same in Fig|4j and the corresponding dash 
lines represent M = p a ■ 47ri? 3 /3. 



where m(r) = J Q p ■ Aitr^dr' . In the above discussions, 
we have got the equations of state, from which we can 
make a further calculation of the mass-radius and mass- 
central density (rest-mass density) relations for quark- 
cluster stars. The results are shown in Fig(SJ for the 
four groups of parameters A to D, and we can see that 
the maximum masses are higher than three times the 
solar mass Mq for all the selected groups of parameters. 
As a comparison, we also plot the mass-radius curves 
for homogeneous spheres with the same central density 
corresponding to each of the four cases. This shows that 
the gravity cannot be negligible only when the stars is 
near the maximum mass, which could be the result of 
the strong self-bounding of quark-cluster stars. 

Conventional quark matter is characterized by soft 
equation of state, and the emerge of quark matter in- 
side compact stars is usually thought to be a reason for 
lowering their maximum mass. The quark-cluster mat- 
ter, however, could have stiff equation of state due to the 
strong coupling. Although the corresponding-state ap- 
proach is just a phenomcnological and empirical method, 
we could still apply it to study the state of quark-cluster 
matter and then understand the observations of pulsar- 
like compact stars. The observed high-mass pulsar PSR 
J1614-2230 with mass 1.97±O.O4M [lj] has received a 
lot of attention, and we can see that the quark-cluster 
stars in our present model could be consistent with this 
observation. Moreover, our model of quark-cluster stars 
could not be ruled out even if the mass of the pulsar 
J1748-2021B (2.74M Q ) in a galactic cluster is confirmed 
in the future. 
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FIG. 7. Data points: experimental data of the reduced melt- 
ing heat H* — H/e. Solid line: fitted curve of eq.(l20|). 

IV. MELTING HEAT 

If the kinetic energy of quark clusters is much lower 
than the inter-cluster potential energy, they may form 
a solid state. We will estimate the latent heat of phase 
transition of quark-cluster stars from liquid to solid state 
by the corresponding-state approach. 

We calculate the ratio of melting heat per particle H 
and e for some ordinary substances [i~9| . and find that 
there is also a good relation between H* = H/e and A*, 
as shown in FigJTl The fitted formula for H* and A* is 

H* = i.i 8e -((A*-o.i2)/i.60) 2 _ (2Q) 

For given N q , e and <r, we can determine A* first and 
then get the value of H* from eq. (j2T))) , thus the melting 
heat H = H*e can be derived. 
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FIG. 8. e and a which determine the same melting heat of 
each cluster. H is chosen to be 1 (black line), 10 (cyan line) 
and 100 MeV (blue line), in two cases N q — 6 (solid line) and 
N q = 18 (dashed line). 



In FigJS] pairs of e and a which determine the same 
melting heat are plotted, where values of H are chosen 
to be 1, 10 and 100 MeV, in two cases N q = 6 and 
N q = 18. The solidification of quark-cluster stars has 
been sug gested to be relevant to the plateau of 7-ray 
burst [20|. and it is found that if the energy released by 
each quark-cluster in the liquid to solid phase transition 
is larger than 1 MeV, the total released energy could pro- 
duce the plateau. We can see that under a wide range of 
parameters in our model, the latent heat could be suffi- 
cient for this way of understanding the plateau of 7-ray 
burst. 



V. CONCLUSIONS AND DISCUSSIONS 

In cold quark matter at realistic baryon densities of 
pulsar-like compact stars, the interaction between quarks 
would be so strong that they could condensate in po- 
sition space, forming quark-clusters, and the stars are 
then called quark-cluster stars if the dominant compo- 
nent inside is quark-clusters. We propose that the in- 
teraction between quark-clusters is like that between in- 
ert gas atoms described by the Lennard- Jones potential, 
and apply the corresponding-state approach to derive the 
equation of state. As a phenomenological and empiri- 
cal method, the corresponding-state approach can avoid 



detailed assumptions of quark-cluster matter as well as 
computation of the many-body effects, and we only need 
to concern about differences between substances. Along 
with of these advantages, there are large uncertainty in 
our results, coming from the Lennard- Jones approxima- 
tion and lack of experimental data source. Even so, the 
corresponding-state approach could give us qualitative 
information about the properties of quark-cluster matter, 
while the exact approach under QCD calculations seems 
to be impossible due to the significant non-perturbative 
effect. Summarily, our two-parameter (e and a) empiri- 
cal approach make it possible to establish a model which 
could be tested by observations. 

The equation of state we have derived by the 
corresponding-state approach could be stiff enough to 
make a star stable even if its mass is higher than 2M@, 
under reasonable parameters. This result is consistent 
with the recent observation of a high-mass pulsar, thus 
the emergence of such kind of exotic matter, "quark- 
cluster matter" , could not be ruled out. The observations 
of pulsars with higher mass, e.g. > 3M Q , would even be 
a support to our quark-cluster star model, and give fur- 
ther constraints to the parameters. Moreover, the latent 
heat released by the solidification of newly born quark- 
cluster stars could help us to understand the formation 
of the plateau of 7-ray burst. 

Whether quark-cluster matter could exist at supra- 
nuclear densities, and what quark-clusters are composed 
of, as well as how to describe their interaction are still 
open questions. On the other hand, the nature of pulsar- 
like compact stars is still uncertain. These are all essen- 
tially related to the non-perturbative QCD, and we hope 
that future astrophysical observations, complementary to 
the territorial experiments, could give us hints to all of 
these questions. 
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